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Abstract 

We study periodic Brownian paths, wrapped around the surface of a cylinder. One 
characteristic of such a path is its width square, w'^, defined as its variance. Though 
the average of w"^ over all possible paths is well known, its full distribution function was 
investigated only recently. Generalising w'^ to w^^\ defined as the iV-th power of the 
magnitude of the deviations of the path from its mean, we show that the distribution 
functions of these also scale and obtain the asymptotic behaviour for both large and 
small w^^\ 



1 Introduction 



Of all the quantities which characterise a near-planar surface, the width, w, defined as the 
standard deviation of the position of the surface, is probably the most frequently studied. If 
this surface is a member of a Gaussian ensemble, then (w^), the ensemble average of w'^, is a 
well-known function of the parameters of the Gaussian as well as the size of the surface. In 
particular, if the surface is a one dimensional object of length T, then (w^) oc T, as T — oo. In 
a recent paper it was pointed out that the entire distribution of w"^, rather than merely the 
average, would provide a better picture of the surface. Indeed, this distribution was shown 
to obey scaling and the universal scaling function was computed. An analogue in critical 
phenomena is the difference between measuring a single critical exponent and finding data 
collapse in a range of parameter space. It is indisputable that the latter gives far more details 
of the system than the former. 

In a similar vein, one may extract further information from the fluctuating surfaces by 
analysing moments or cumulants other than the second. For Gaussian surfaces, the averages 
of these can be related to (w^), but their distributions are less trivial. Unlike the case for w'^, 
the Laplace transforms of these distributions are associated with non-Gaussian field theories. 
Equivalently, for one-dimensional surfaces, we must now deal with quantum anharmonic oscil- 
lators. In this paper we generalise the analysis begun in and show that these distributions 
also scale, while their universal functions can be found in the asymptotic regimes where the 
argument is small or is large. 

In order to have sufficient analytical control over the fluctuations, the one-dimensional 
surfaces in [Q were modelled by Brownian paths {ht} in the "time" interval < t < T, with 
periodic boundary conditions. Thus, a surface may be thought of as a line wrapped around a 
cylinder of circumference T. By identifying the continuous path ht as the height of a surface, 
the "width square" is a simple functional quadratic in hf. 



w'^[ht\ = ipl = ht-ht 
where the average, /, of a function / is defined as 

7{ht) = ^l^dtf{ht), (2) 

and where ipt = ht — ht is the usual deviation from the "average height". Being Brownian 
paths, the weight of each ht is a simple Gaussian, i.e., exp {— Jq^ dt ^h^ } , where ht = dht/dt. 
For the ensemble average, denoted by (...), one must sum over all appropriate paths, i.e., 
perform functional integrals (or path-integrals 0) over periodic functions ht. Since none 
of the quantities we encounter will depend on the overall height, we may just as well sum 
over {ipt}, with special attention paid to the constraint dtipt = 0. With this set up, it is 
standard to show that {w"^) = T/12. 

The key to calculating the probability density P{w'^\ T) was expressing it the form of a 
path- integral and noting that its Laplace transform can be computed explicitly [jll. Beginning 
with 



P{w'-T) = {6[w^-\ht~ht'' 



U jv[h]5{w^-[ht-ht]^Yyi^{-^h1^ , (3) 
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where A/" is a normalisation constant, one can easily show that it has the scaling form 
P{w'^;T) = {w'^)~^^{x), where x is the scaling variable w'^/{w'^). Inverting the exact re- 
sult for its Laplace transform 



/AT/2 

G{X;T)= ' ' , (4) 
sinh(JAT/2) 



the associated scaling function is obtained: 



$(^) = yE(-ir"'^'exp --n^a; , (5) 



n=l 



from which the asymptotic properties can be found. 

The aim of this paper is to extend these results to arbitrary absolute central moments: 



w 



i^)[h,]^\i,^\^ = \ht-ht\^ . (6) 



Notice that, for even N , these are precisely the central moments associated with a particular 
configuration ht. Using the absolute values of ip, it is possible to study arbitrary > 0, and 
these are known as the absolute central moments (associated with a specific path). The outline 
of the paper is as follows. In section 2, a brief, self-contained, summary of the formalism for 
this general case is included. Then, we show that Pn{w^^^;T) does indeed scale. Analytic 
expressions for the asymptotic behavior of the scaling function are found, for large and small 
arguments, in sections 3 and 4 respectively. We conclude with some general comments in 
section 5. 



2 Scaling form for the probability distribution 

Following Ref. |^ , let us express the full distribution of w^^^ as 

PNiw^""^; T)=M j V[h] 6 - Jfh^W" ) exp 
Its Laplace transform, which generates the moments of (^, is 



T ^ 

-hi 
2 * 



(7) 



dCP^{C-T)e 



M V[h] exp 



'AC 



Tt 



hi -\\ht~htW 



Writing this expression in terms of ipt and using the condition G'Ar(0; T) = 1 to determine the 
normalisation constant M , we obtain 





exp 


-lodt 






JV[^]6{j^dtij,) 


exp — Jq dt 







(9) 



Apart from the constraint dtipt = 0, this is the imaginary time propagator for a quantum- 
mechanical particle moving in the potential IV'I^- As we will see in section 4, this connection 
will be exploited in the study of the asymptotics of Pn for small w^^\ 
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To elucidate the general structure of Pn and its scaling properties, it is useful to define a 
dimensionless variable t = t/T and then to make the following rescaling: 



Defining 

we have, for the generating function 



IV[x]6{j,'dTXr) 


exp 


- Jo dr 


ly2 I A 
2At ^ N 


Xr 


n' 


JV[x]6{j^dTXr) 


exp — Jq dr 


2 At 





(10) 



'11^ 



(12) 



where Xt now means dxr/dr. This expression shows that Gn{^', T) is a function of only, so 
that we may write 

f;^(/i) = G;v(A;T) (13) 
Expressing Pn(w^^^;T) as the inverse Laplace transform of g: 



PNiw^^'^^T) = r ^^^.(iVAT^/2)exp(A^(^)) 

^ dfi 



-too 

1 



]S[rpN/2 



27Ti 



qn^Ij) exp 



/i w 



(14) 



we see that T^/^P, 



(N). 



we simply note that (w^^^) is proportional to T^/^: 



T) is a function of w^^^T^^"^ only. To put this in the scaling form, 



dCCP^iCT) 

dGi\i 



dX 



A=0 
d 



djjL 



(15) 



Ai=0 



Using this result, (|T^ may be written as 

P^(^(^);T) 

where 



and 



— ^^(/i) exp(/ia;|c/^(0)|) 

-ioo ZTTi 



(16) 
(17) 



Equations (p^ -(|T7D are an expression of the fact that Pn scales. From (pTSf ) one sees that 
in order to find an explicit form for $Ar(a;), the function gN{iJ) has to be determined. For 
N ^2 this can only be determined in special limits. Nevertheless, as we shall see in the next 
two sections, one can go quite a long way in finding the limiting forms of $Ar for large and 
small X using path-integral techniques. 
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We end this section by noting that, since we have Brownian paths, the ensemble averages 
of the moments (w^^^) are quite simple. For even N, they are (iV - l)!!(r/12)^/2. For 
arbitrary A^, we rely on standard analytic continuation techniques to arrive at 

|,;(0)| = r(^)/{6^/W/^}. (19) 

However, we will continue to write 5'^(0), both for convenience and to emphasise its role. 



3 The scaling function for large x 

It is not surprising that the behaviour of $iv(x) for x 1 is controlled by the singularities 
of 5'Ar(/i) for near the origin of /i. In the N = 2 case, these singularities are simple poles on 
the negative real axis. In general, we expect a branch cut and its discontinuity near fi = 
to dominate the asymptotics of $7v One approach is to compute these singularities, using 
techniques of asymptotic expansion around instanton solutions 0, followed by an inversion 
of the Laplace transform. Alternatively, one can insert (0) into (|T8|) , find the saddle point 
in the combined space of {/x, Xt} and perform the Gaussian integrations over small variations 
in its neighbourhood. We shall follow the latter route, which seems to be simpler. 

As will be seen, we are able to compute explicitly only the leading two terms in In ^j^, 
i.e., x^/^ and In x. Thus, we will drop all proportionality constants, for clarity, when writing 
$Ar. In particular, the denominator of (0) is clearly x independent and will be suppressed, 
though its role in regulating the Gaussian integrations is obviously essential. So, consider 

$Ar(x) oc J dfiV[x] 6 (^J^ drxr^ exp l^fix\g'^{0) \ - dr 

The attentive reader may object to the factor i in ( pISD being dropped as well, since $Ar should 
be real. A careful analysis will reveal that the functional determinant carries this factor, i.e., 
the contour for fi near the saddle point will be parallel to the imaginary axis. 

From (|20|) , we first seek the saddle point, to be denoted by {— M, X,-} in {/x, Xr} space. 
Note that we have anticipated this point to be on the negative real /i axis, so that M will 
turn out to be a real and positive. Setting to zero the variation of the exponent in (^) with 
respect to /i and Xr , we have the equations for the paths which dominate the path-integral: 

Nx\g'j,m= ['dT\X,\^, (21) 
Jo 

and 

Xr + Msgn{Xr) |X,|(^-^) = . (22) 

Note that, in order to ignore the singularity at X = 0, we must restrict our attention to X > 1 
here. The latter being an equation of motion for a classical particle of unit mass moving in a 
potential V{X) = f |X|^, these can be solved easily. The importance of M > is also now 
evident, since periodic Xr's would have been otherwise impossible. Finally, the constraint 
lo drXr = can be satisfied trivially. 
Integrating (|2^) once, we have 

l^r + ^\Xrr = e, (23) 



iy2 I R\y \N 



(20) 
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where e is a constant, representing the total energy in the mechanical analogy. Assigning 
r = to the point of maximum amplitude (and so, zero velocity), we set 

^=^^r (24) 

where Xq > 0. Equation ( ^31) can be integrated again in the usual manner, exploiting both 
the periodic boundary conditions and the constraint. Clearly, we need to focus on only a 
quarter of the period, so that 

1 ^ x('-^)/'(iV/2M)i/2 di{l - er^'^ , (25) 
4 Jo 

where C,t = Xt-/Xq. The integral is proportional to the beta function B{j^, |), which we will 
denote simply as B. Now we have a relation between M and Xq: 

Inserting this into Eqn. (|2T|) , eliminating dr in favour of dX/X and using (^31) , we have 

Nx\g'^m = 4/ {dX/X)\X\^ (27) 
Jo 

Using B{1 + a, P) = B{a, P)-^^, we find both Xq and M in terms of the scaling variable x : 

N 



X, 

and 



o = {iV(y + l)a:|^^v(0)|} (29) 



M= (^^j|iV(- + l)a;|<(0)|| . (30) 

With all parameters of the saddle point explicitly determined, we evaluate the exponential 
in ([20|), i.e., the total "action". The result is —MNx\g']y{0)\/2, so that the leading asymptotic 
behaviour of $Ar is 

<l>^(x)~exp{-5^x2/^} , (31) 

S^ = AB'[NC-^ + l)j k5v(0)|(^/^^ (32) 

For N = 2, B = 71 and, from (|l^), |fi'Ar(0)l = 1/24, so that we recover the result in (|). 

Next, we turn to the computation of the prefactor, for which it is necessary to study the 
Gaussian fluctuations about the saddle point (— M, X^). 

To accomplish this, we write /i = — M + fi and Xt = X^ + Xt in (|20D and shift the 
integration to (/i, x^). Expecting the integrals to be dominated by small (/i, Xr), we expand 
the argument of the exponential to second order in these fluctuations. The zeroth order is 
given above in (^) while the first order vanishes by the choice of the saddle point. At the 
second order, the result is the quadratic form 

-fi f'dTSgniX,)\Xr\^^-''^Xr + l f'drx^MXr, (33) 

JO 2 Jo 
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where M. is the operator 

dr 



M = -—-M{N -l)\X^\'^''-^K (34) 



Note that, strictly speaking, we should impose > 2, so that the singularity at X = can 
be ignored here also. 

Before carrying out the integration over (/i, Xr), we discuss several important points. 

Firstly, note that is a hermitian, Schrodinger-like operator. It can be diagonalised, and, 
since our problem is based on a finite interval with periodic boundary conditions, it has a real, 
discrete spectrum with real, periodic eigenmodes. Thus, x can be expanded in terms of these 
modes and the the functional integral can be defined as over the amplitudes in this 

expansion. We will demonstrate that there is a single zero eigenvalue, associated with a zero- 
mode, which must be handled with some care. Defining to be the operator ( p4|) restricted 
to the subspace orthogonal to this mode, we will argue that the spectrum of A4' is positive, so 
that its inverse is well defined and the Gaussian integration in (|33|) over these modes is simple. 
Furthermore, the product M|X^|(^-2) in (||) can be written as {MX^^~^'^)\^r\'-^-'^\ But, 
according to (^), this quantity is independent of x, so that the spectrum of both M. and A^' 
are also x- independent. Thus, the result of the Gaussian integration, which involves detA^', 
is also x-independent. Since we are dropping all proportionality constants in this study, this 
factor will be neglected. 

Secondly, we will show that the zero- mode is also absent from the off-diagonal fi - x part. 
As a result, at the quadratic level, the x integration leads to a distribution for fi of the form 
exp (P/i^/2), where 

V = I drdr' |X,|(^-i)Al'-i|X.,|(^-i) (35) 

is positive. To make sense of performing such an integral (over /t), recall that the contour in 
the /t-plane is parallel to the imaginary axis, by definition of the inverse Laplace transform. 
We have simply chosen to have it run through a saddle point on the real axis: yU = — M. Thus, 
it is entirely consistent to choose fi to be pure imaginary, so that not only is the Gaussian 
integral well defined, but also the right hand side of (|20D is now real. Further, we may again 
use Xt = Xq^t- to extract the sole dependence of V on x, given that neither Ai' nor are 
functions of x. As a result, the integration over fi leads to a factor proportional to 

1/VV oc X^'-"^) . (36) 

Thirdly, for Schrodinger problems in one dimension, the energy levels are non-degenerate, 
unless its potential is a constant. The implication for us is that, except for the N = 2 case, the 
spectrum of Al is non-degenerate. Further, since our interval is finite, the spectrum is discrete, 
with eigenvalues increasing monotonically with the number of nodes in the eigenfunction. 
Now, the constraint 6 Jq dr Xr) must clearly be satisfied by both x and the eigenfunctions. 
Therefore, the latter must have an even number {2n > 2) of nodes. We will show that the 
zero-mode has the lowest allowed n. Since it is not degenerate, the rest of the spectrum of 
Al, i.e., the spectrum of Ai', must be positive. 

Now, let us provide some details of the zero-mode. Its origin is the translational invariance 
in our problem. To identify it explicitly, differentiate (^) once and find 

M^ = 0. (37) 
dr 
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Thus, Xr is an eigenfunction of Ai with zero eigenvalue, so that it is clearly the zero-mode. 
Next, to show that it is also absent from the first term in (|33|) , we note that sgn(X^)|X,-|'^^~^'' 
is actually dV{X)/dX. Therefore, if Xt is then the integrand is a perfect derivative of a 
periodic function and this term vanishes. The conclusion is that the zero-mode is completely 
absent from the quadratic from (pSl). Thus, for the Gaussian approximation to be valid this 
mode has to be treated separately. The standard technique to deal with this situation is the 
method of collective coordinates [Q. The symmetry here is translational invariance of ( pi]) 
and (^2]) , which implies that the particle (in the mechanical analogy) can be at the point Xq 
at any other time tq G [0, 1]. Thus, there is a one-parameter family of solutions, labelled by 
To, all of which satisfy the differential equation (p^) and the constraint Jq drXr = 0. Indeed, 



the zero- mode X^- can be recognized as the difference between X^. and X^-.^q , to lowest order 
in Tq. Using this parameter instead of the amplitude of X^- in the normal mode expansion 
of X enables us to sum over this component of x, even though it is entirely absent from the 
weights. 

Having made these remarks, we may now evaluate the Gaussian integral over x. In 
principle, we would expand x terms of all the eigenfunctions of Ai, except the zero- mode 
(associated with the index n = below). Therefore, Xr = Sn^o '^"Xi"'' • the collective 
coordinates method |Q, the integral over the functions x is replaced by integrals over tq (which 
"replaces" 0) and {ctnl?^ 7^ 0}. The Jacobian of this transformation is simply which 
is proportional to Xq and carries the sole x dependence here. Now, the integration tq is 
trivially unity. Meanwhile the integration over the rest of the modes ({a„|n 7^ 0}) has already 
been discussed, in connection with the first remark above. No x dependence appears here. 
Finally, the integration over fi leads to (|36D. Summarising, the Gaussian integrals yield a 
prefactor for (pT]) proportional to Xq^ ^\ Using (pOj), we have our final result: 



TV — 2 

ln$^(s) ~ -^Tva:'/^ ^ln(x) + 0(l) (x > 1) , (38) 



where, explicitly. 



''^-3iv(iv + 2)l,r(^)j { j ■ 

Although there are some technical limitations to our derivation for for the N = 2 case, 
they can be circumvented, so that this result is in fact valid for N > 2. 



4 The scaling function for small x 

To arrive at $iv in the other limit: x -C 1, we will use a very different approach, which we will 
first sketch briefly. As in the a; ^ 1 case, we are able to compute only the two leading terms 
in In $ so that all proportionality constants will be dropped. Unlike in the previous section, 
we will first compute the asymptotic behaviour of Gn{X;T) for large A and then rely on a 
saddle point method to find $Ar. Our starting point is the unsealed version of Gn{X] T) given 
by (^. So first consider the numerator of the right-hand-side, but without the constraint 
dtipt = 0, and define 

a = \/T (40) 
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Then, 



exp 



dt 



\N 



OC 



(41) 



where is the (n+l)-th eigenvalue of the evolution operator which consists of the quantum- 
mechanical Hamiltonian for a particle moving in the potential As T ^ oo, the ground 
state dominates, leading to the Feynman-Kac result that the right-hand-side of (|4l|) can 
be replaced by e~^o with only an exponentially small error. Now, the dependence of 
on a can be extracted through dimensional analysis alone, since a has dimensions of inverse 
time to the power ^(A^ -|- 2), while Eq^"^ itself has dimensions of inverse time. Thus, we have 

Eq^^ oc a^/*^^"'"^'' and arrive at the the leading behaviour of the path integral in (|4l|) , for 

T > a-2/(iv+2)_ 

Having outlined the essential idea, let us now be more systematic. First, to implement the 
constraint Jq dtipt = 0, we write the usual integral representation for a delta function. Next, 
instead of normalising the path integral as in (^, we use (j'2(A;T) as the denominator. This 
avoids the complications from the a = system, which lacks a discrete spectrum for a simple 
application of the Feynman-Kac result. On the other hand, 6*2 (-^i T) is known explicitly (see 
d^)), so that the ratio Gn/G2 is essentially our goal. Thus, we consider, from (|), 



G2{X;T) 



JZodoo JV[ij]exp-j;[dt 


+ a\^t\^ + iujiJt 


JZ.duJV[^p] exp-/5^rft 







(42) 



Of course, the presence of iojip is rather unusual. However, as we will see, its sole effect is 
to give a non-trivial prefactor to G in this limit. We first focus on the denominator, which 
consists of simple Gaussian integrals. Without the constraint, the iuil) term is absent, so that 
we may use results from the quantum mechanics of a simple harmonic oscillator to obtain 

G2 ~ e~v^^"^. Comparing with (1), we see that this is indeed the leading behaviour. The 
next leading term, contained in the prefactor, would arise from the constraint. In this case, 
it is easily dealt with, by defining ip[ = ipt + iuj/2a and factorizing the integrand. The ip' 
integral gives us the previous asymptotic form, while integration over uj produces the desired 
prefactor: ^ajT . For the N ^ 2 case, there is no similar luxury of factorization, of course. 
But, by reinterpreting the above steps, i.e., computing the functional integral with the iujil) 
term, as finding an (^-dependent ground state energy: i?Q^'*(a, cu), a way forward can be seen. 
Thus, we write 



j Vli)] exp-^ dt 



oc 



(43) 



with 



E 



(2), 



(«,cu) = ^S'^(«,0)+cuV4a. (44) 

It is now clear that the integral over u gives the prefactor. Another route to the understanding 
of ( ^^ is to consider adding a real, linear potential: hip first. Clearly, the ground state energy 
is well defined, allowing us to arrive at the above result by analytic continuation to pure 
imaginary h. 



Generalizing this approach to the N ^ 2 case, we write the numerator of (^21) as 

-i?r'(a,a;)T ^ 



dtu } 



^^^-Er(a,.)T 



(45) 
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Unlike the = 2 result (||), eJ^^ {a,uj) does not have a dependence on lu which terminates 
at second order. It is clear, however, by expanding the exponentials in the analogue of (E3|) 
for general N, that it has the general form 



SN) 



(46) 



So since cu has dimensions of time to the power —3/2, e^\a) has dimensions of time to the 
power (3m — 1). Therefore, 



,W(c,)=e^W(l)a-2(3™-i)/{^+2). 



(47) 



To see that only the first two terms in this expansion are relevant, recall that the inte- 
grand depends on EQ^\a,uj)T only. Changing the variable of integration from u to oo' = 
T^/^a~^/('^+^)a;, so that the quadratic term in (|46|) is independent of a and T, we see that the 
m > 2 terms are of order l/[a^/^^"^^''T]''"~^. Since this is proportional to l/^2(m-i)/(Ar+2) 
we are interested in the yU ^ 1 limit, we are justified in neglecting all terms beyond m = 1. 
Focusing on the first two terms in (^6|), we write 



Ef)(«,u;)T = £:a2/(^+2)2^ ^ ^/^^.^2 ^ 



(4J 



Note that S = eQ^\l) is just the zero point energy of a particle in the potential and is, 
therefore, positive. Meanwhile, had oo been pure imaginary, E^^-' would certainly have been 
lowered, implying that S' is also positive. We note, parenthetically, that it can be estimated 
in our approach, since it is just the correction to the ground state energy in second order 
perturbation theory. Having S' > 0, the integral over u' is well defined and provides a simple 
constant. 

Inserting (HSI) into (53), we have 



Using this, together with a similar expression for the denominator, 

" ^^g-£;,f)(a,.)T ^ ^1/2^-1/2 \a'/^T/V2 



_^^2/(Jv+2)y] |i + o(a-2/(iv+2)2^)| 



(49) 



(50) 



we find 



Gn{X; T) A2/(^+2) t-2/(^+2) exp 


'_gy2/{N+2)rj.N/{N+2) 


G2(A;T) AV2 r-V^ ^^p 


-AV2Ti/2/y2' 





{l + 0(AT^/2)- 



-2/(Af+2)- 



where we replaced the a's by the original A/T's. Finally, using the exact expression 
G2, we arrive at 

G;v(A;T) OC A2/(^+2)y^/(iV+2) ^^p J_^;^2/(iV+2)y^/(iV+2)] { 1 + Q ( AT^/^) ^2/(^+2) | ^ 

Alternatively, we can write it in terms of the scaling form gN^fJ^) introduced in section 2: 



(51) 
) for 

(52) 



g^ifl) OC /i2/(^+2)^^pr_^(^/^)2/(7V+2)l|^^Q^^-2/{iV+2))|_ 



(53) 
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This is the form which Gn or takes in the /i ^ 1 hmit. It should be contrasted with the 
approach used in the last section, which was justified in the regime where <C 1. Thus the 
results of the last section are complementary to those derived here. 

All that remains to find is to take the inverse Laplace transform of (0) : 



oc 



djj, jj, 



2/(N+2) 



exp 



Defining u = x^^'^'^^^^ fi gives, for the exponent in (0), 



(54) 



(55) 



so that the integral may be evaluated by steepest descent if x <^ 1. Since S > 0, there is an 
extremum at a real positive value of u, which is of order unity. Thus, the corresponding value 
of /i is of order a;-(^+2)/A^ ^ justifying all the approximations we made above, based on 
^ oo. Carrying through the integration gives the asymptotic behaviour: 



In $jv(x) ~ -JCnx-^/"" - ln(a;) + 0(1) {x < 1) 



(56) 



where /Cat is a positive constant which depends on some explicitly known functions of and 
on S, the ground state energy of a quantum mechanical particle of unit mass in a potential 
^(C) = ICI'^- Unfortunately, the last item is not known for general N, so that we must be 
content with just an implicit expression. Nevertheless, conjecturing that it is monotonic in 
A^, we can place a bound on it: vr^/S, which is the result for A^ — > oo, where the particle is 
confined to |C| < 1 by an infinite square well. Since £ is I/V2 for A^ = 2, it is fairly well 
bounded numerically even though its precise value is unknown. Finally, we note that (|56D 
reduces to the known result given in ref. []T| for the N = 2 case. 



5 Conclusions 

To summarise, we have generalised the study of width distributions of a Gaussian path (ran- 
dom walk) to include moments higher than 2. Unlike the N=2 case, we are unable to find 
closed form expressions for the Laplace transform of these distributions. Nevertheless, we 
can compute their asymptotic behaviour, for both large and small arguments: eqns. ( pH] ) and 
(|56|). Here, we conclude with a few remarks. 

Firstly, note that, though the averages of the even N > 2 moments are trivially related to 
{w"^), their distributions are not so simple. In particular, we believe that there is no way to 
obtain the general $Ar or from the N=2 result. They contain different information about 
the paths. For example, the following two paths lead to the same w but distinct u;W: (i) 
^/'t = 1 for < t < T/2, -1 for T/2 <t <T and elsewhere, vs. (ii) V'f = 2 for < t < T/8, 
—2 for 7T/8 < t < T and elsewhere. Thus, only paths with short excursions to large ip^s, 
which are presumably rarer, contribute to the large x tail of $Ar with large A^. Presumably, 



this is reflected by eqn. (|31|). 

Secondly, we have seen, from Feynman's path-integral formulation of quantum mechanics, 
that the Laplace transform of the distributions for the A^-th absolute central moments of 
random walks are intimately related to the propagation kernels of a quantum mechanical 
particle confined to a potential of the form V{() = \C\'^ ■ As a result, the behaviour of ^n{x) 
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in the large x limit should be controlled by the properties of Gat (A) for small A. We note that 
the paths which dominate Gn here are the ones near the "classical path" (^). Meanwhile, 
for the opposite limit, small x and large A, the dominant contribution to Gn comes from the 
ground state, which is, in a sense, "the furthest from classical". It is interesting that these 
opposing aspects of the quantum mechanical problem find their way into the opposite ends 
of the asymptotics. 

Thirdly, let us recall that these distribution functions are supposedly universal, in the lan- 
guage of renormalisation group. In other words, they should capture the large-scale properties 
of a one-dimensional object stabilised by non-zero tension, independent of the microscopies 
of the system. Thus, we may expect the same behaviour from an interface between different 
phases, a polymer in ci = 2, or a simple model of random walk. The last of these is particu- 
larly easy to study, using Monte Carlo simulations. In this connection, we should caution the 
reader on the N oo limit. Simulations necessarily deal with systems with finite L. To ob- 
serve universal properties, we must let L approach infinity. However, this does not commute 
with — >■ CX3, so that the results above should be used with some care, if large values of 
were to be used in the investigation. 

Finally, it is natural to speculate on further generalizations of this study. Obvious candi- 
dates are interfaces in higher dimensions and those controlled mainly by curvature terms 
H^. Both will involve single component fields. On the other hand, we could investigate ran- 
dom walks imbedded in higher dimensions, such as physical polymers in 3-d solutions, where 
the periodic boundary condition would correspond to ring polymers. In this case, we would 
need multi-component fields and an appropriate generalisation of the concept of "width" . In 
the simplest scenario, P(w^) will be nothing more than products of the distribution in 
However, it is clear that P/v will be much more interesting! While the formulation of any of 
these problems is trivial, we anticipate that the analysis itself will be be rather difficult. 
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